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We study static and dynamical critical phenomena of chiral symmetry breaking in a two-flavor 
Nambu-Jona-Lasinio constituent quark model. We obtain the low-energy effective action for scalar 
and pseudoscalar degrees of freedom to lowest order in quark loops and to quadratic order in the 
meson fluctuations around the mean field. The static limit of critical phenomena is shown to be 
described by a Ginzburg-Landau effective action including spatial gradients. Hence static critical 
phenomena is described by the universality class of the 0(4) Heisenberg ferromagnet. Dynamical 
critical phenomena is studied by obtaining the equations of motion for pion fluctuations. We flnd 
^— ■») . that for T <Tc the are stable long- wavelength pion excitations with dispersion relation lutt (k) — k de- 

' scribed by isolated pion poles. The residue of the pion pole vanishes near T^. as Z oc l/[ ln(l — T/Tc)] 

and long- wavelength fluctuations are damped out by Landau damping on a time scale trei(fe) oc 1/fc, 
Sh , reflecting critical slowing down of pion fluctuations near the critical point. At the critical point, the 

, ^ ' pion propagator features mass shell logarithmic divergences which we conjecture to be the harbinger 

of a (large) dynamical anomalous dimension. We flnd that while the classical spinodal line coincides 
, with that of the Ginzburg-Landau theory, the growth rate of long-wavelength spinodal fluctuations 

• has a richer wavelength dependence as a consequence of Landau damping. We argue that Landau 

damping prevents a local low energy effective action in terms of a derivative expansion in real time 
at least at the order studied. 
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■ I. INTRODUCTION 

O ; 

i-S^ ' The spontaneous breaking of chiral symmetry in the QCD vacuum has fundamental significance in understanding 
^H, the nonperturbative nature of hadron dynamics. In QCD with two massless quark flavors at zero baryon density, 
O I the chiral phase transition is expected to be second order and to be described by the universality class of the three- 
D , dimensional 0(4) Heisenberg model 1', Recent lattice simulations for two- flavor QCD in the chiral limit at finite 

■ temperature and vanishing baryon chemical potential seems to indicate that chiral symmetry is restored above a 
J> I transition temperature To ~ 170 MeV, whereas the demonstration of the expected 0(4) universality class still is 

somewhat ambiguous 0]. 

In ultrarelativistic heavy ion collisions the current theoretical understanding suggests that a thermalized and almost 
^ 1 baryon-free quark-gluon plasma (QGP) may be formed at a time scale of order 1 fm/c with an initial temperature 
■ " ' larger than the transition temperature. This quark-gluon plasma then expands hydrodynamically and cools almost 
adiabatically, until the transition temperature is reached at a time scale 10 — 50 fm/c depending on the initial 
temperature Q- At BNL Relativistic Heavy Ion Collider (RHIC) and CERN Large Hadron CoUider (LHC) initial 
temperatures of the quark-gluon plasma as high as T ~ 300 and 500 MeV, respectively, are expected to be achieved in 
the most central collisions, hence providing the experimental possibility to study the restoration of chiral symmetry 
and its subsequent spontaneous breaking as the plasma expands and cools. 

While lattice QCD furnishes a nonperturbative tool to study the thermodynamic equilibrium aspects of the chiral 
phase transition, the nonequilibrium aspects which are of fundamental importance and ubiquitous in the formation 
and evolution of the QGP in ultrarelativistic heavy ion collisions cannot be accessed with this approach. They have to 
be studied within the framework of nonequilibrium field theory. However, due to the rich nonperturbative structure 
of QCD current theoretical investigations are limited only to models that incorporate the relevant symmetries and 
low-energy effective field theories of QCD. Some popular models of chiral dynamics are the linear sigma model 
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the Gross- Neveu model '6| , and the Nambu- Jona-Lasinio model 'l] . The thermodynamic equilibrium aspects of these 
models at finite temperature and/or quark chemical potential have been studied extensively in the literature [sLI^ITol 
HH l , in particular lattice gauge theory has provided an impressive body of results on the phase diagram for temperature 
and chemical potential (see for example ITsl Il5l| . Nevertheless, a microscopic understanding of their dynamical 
nonequilibrium aspects has just begun to emerge. 

It has been recently argued that nonequilibrium dynamical phenomena near the critical point such as spinodal de- 
composition, critical slowing down and formation of correlated pion domains "31,11611 may lead to distinct event- by- event 
observables in the pion distribution ^17.j as well as in the low-energy spectrum of photons [T^ . The nonequilibrium 
dynamics in the linear sigma model has been used to study the possible formation of large correlated domains of 
disoriented chiral condensate (DCC) in ultrarelativistic heavy ion collisions |l6l |. Thus non-equilibrium phenomena 
associated with chiral symmetry breaking may have important observational consequences in ultrarelativistic heavy 
ion collisions. 

While several studies of the nonequilibrium dynamics focused on either a discrete chiral symmetry HOl or a stochastic 
description of long- wavelength fluctuations above the critical point in the chirally symmetric phase |20| , the dynamics 
of long-wavelength fluctuations at or below the critical point has begun to receive attention only recently. 

Critical slowing down and the relaxation of long-wavelength fluctuations at or near the critical temperature has 
been studied in the large- iV limit of the 0{N) linear sigma model as well as by implementing a finite-temperature 
renormalization group in Ref. |i2l| . The results of the rcnormalization group indicate a new dynamical critical exponent 
z ^ 1 + 1/27 which results in a dispersion relation for long- wavelength critical fluctuations of the form u!{k) oc 
with vanishing group velocity and critical slowing down in the long- wavelength limit |2lj| . This study focused on 
long-wavelength fluctuations for T > but did not address the relaxation dynamics of pions below the critical 
temperature. 

In Ref. the real-time propagation of pions near the QCD chiral phase transition was studied in an effective 
chiral theory which in principle is obtained by integrating out the quark (and gluon) degrees of freedom of QCD. The 
effective low-energy theory proposed in Ref. 22| on the basis of symmetries, is local in the sense that the terms that 
describe the space-time variations of the pion fleld are written in a gradient expansion up to second order in space 
and time derivatives. Based on this effective chiral action, the conclusion presented in Ref. |22j was that the real-time 
propagation of soft pions below the critical point can be solely expressed in terms of static correlation functions. An 
important result of such analysis is that the pion velocity decreases near the phase transition. While the results of 
this study apply to temperatures below the critical and cannot be directly compared to those of Ref. (which 
apply either at or above the critical temperature), there does seem to be a discrepancy between these results at least 
if the results of Ref. '2^ are extrapolated to the critical point. In particular the study in Ref. |^ suggests that 
there is no local description of the dynamics in the low-energy effective theory, namely the space-time variation of 
long-wavelength fluctuations cannot be described in terms of a low-energy effective theory with a simple derivative 
expansion for the terms that determine spatiotemporal variations. In the analysis of Ref. |2l| this breakdown of locality 
near the critical point emerges in the form of logarithmic singularities near the mass shell of the single (quasi)particle 
excitations, which a resumniation via the renormalization group translates into the anomalous dimension z |2l| . 

Goals of this article. In Refs. 0,1121 the starting point of the study was a scalar field theory with a Lagrangian 
density that describes the spatiotemporal variations in a local manner in terms of a derivative expansion. In Ref. |2l| 
the field theory studied was the linear sigma model, while in Ref. [22l | it was the nonlinear sigma model. In both 
cases these are low energy effective field theories where the terms in the Lagrangian that describe the spatiotemporal 
variations are local derivative terms. 

In this article we investigate the non equilibrium aspects of the chiral phase transition below the critical temperature 
in the chirally broken phase in the two-flavor Nambu- Jona-Lasinio (NJL) constituent quark model in the strict chiral 
limit of massless (up and down) quarks. In this model scalar and pseudoscalar mesons arise as bound states of 
quarks and are not present as elementary fields in the Lagrangian. Their dynamics is completely determined by the 
underlying dynamics of quarks. The effective action for the meson fields is obtained systematically by integrating out 
the quark fields, thus the low-energy effective action for scalar and pseudoscalar mesons can be obtained explicitly 
and systematically. 

Our main goals are the following: 

(i) To obtain the low-energy effective action up to quadratic order in the scalar and pseudoscalar fluctuations 
around the mean field in the critical region. Namely, we extract the Ginzburg-Landau effective theory, both in 
the static limit as well as in real time. This is achieved by integrating out the quark degrees of freedom up to 
one-loop level. In the NJL model with Nf fiavors and Nc colors, the one loop effective action is exact in the 
large NfNc limit j^, |9i . Thus although we will focus our discussion on the phenomcnologically relevant case of 
two fiavors and three colors, the one loop result can be argued to be exact in the limit of large NfNc. 

(ii) To study the real-time dynamics of long-wavelength pion fluctuations in the chirally broken phase near and at 
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the critical temperature. In particular we focus on extracting the dispersion relation for pions as well as their 
relaxation. Furthermore, we also address the dynamics of long-wavelength spinodal fluctuations and extract 
their growth rate in the long- wavelength limit. 

(iii) To establish the validity of the assumption on the universality of the critical phenomena. The SUv(2) x SUa(2) 
symmetry of the NJL model suggests that the low-energy effective theory for scalar and pseudoscalar mesons 
will be in the same universality class as the three-dimensional 0(4) Heisenberg model 0. We study this aspect 
both in the static as well as in the real-time effective action up to the one loop level in quark loops which is 
exact in the large N^Nf limit |^. 

The NJL model is clearly not QCD, to begin with it includes the effect of gluons as an effective four-Fermi 
interaction and confinement is not a property of the theory. These shortcomings of the model notwithstanding, an 
exhaustive study of its properties point to its remarkable successes in describing the phenomenology of the 
low-energy degrees of freedom, namely the pions. In particular important consistency relations for the two- flavor case 
such as the Goldberger-Treiman and Gell-Mann-Oakes-Renner have been shown to be obeyed by this model (for a 
detailed analysis see |^|3)- The coupling (and cutoff) of the model are fixed so that the pion decay constant at zero 
temperature is reproduced The NJL model being solely a (constituent) quark model, the meson fields emerge 

as quark-antiquark bound states. Hence the kinetic terms of the meson fields in the effective low-energy field theory 
are completely determined by quark loops. Thus the static and real-time aspects of fluctuations beyond the mean 
field are completely determined by the correlation functions of quark operators. 

This article is organized as follows. In Sec. ^1 we briefly review the two-flavor NJL constituent quark model and 
to focus our comparative study. For comparison, we review as well the Ginzburg-Landau description of the 0(4) 
linear sigma model. In Sec. IIIII we obtain the meson self-energies to lowest order in the quark-loop expansion and 
establish the low-energy effective action in the static limit. In Sec. lIVI we study the dynamics of long-wavelength pion 
excitations near the critical point by obtaining and solving the equations of motion for small amplitude fluctuations 
around the mean field configuration. Our conclusions, summary of main results and further questions are presented 
in Sec. |V| 



II. THE NAMBU-JONA-LASINIO MODEL AND THE LINEAR SIGMA MODEL 

The constituent quark model that we consider in this article is the NJL model in the strict chiral limit (vanishing 
current quark masses) for Nc colors and Nf flavors. This model is defined by the Lagrangian density 

£njl = y^i^i' + ^[(^V-)' + i^il^r^pY], (2.1) 

where ip is column spinor describing Dirac quark fields with Nc colors and Nf flavors, G is the dimensionful four- 
fermion coupling, and t* (j = 1, 2, 3) are proportional to the Pauli matrices and normalized as 

Ttt't^ ^ Nf6'^, (2.2) 

with Nf the number of quark flavors. This model has been studied in detail 0| and found to describe fairly well 
the low-energy phenomenology of pions, which are the lightest hadronic excitations in QCD. The Lagrangian density 
-Cnjl is invariant under axial-vector SU(2)y x SU(2)^ transformations and can be regarded as the quark part of 
a low-energy effective theory of QCD after integrating out the ghion degrees of freedom. The contact four-fermion 
interaction is not confining and renders the NJL model non-rcnormalizable. A regularization scheme must therefore 
be employed to cutoff the divergent quantities that occur, and the model must be understood as a cutoff effective 
field theory. 

In the literature, various regularization schemes have been adopted in the NJL model 0,1^. In the finite-temperature 
case under study the noncovariant and physically intuitive three-momentum cutoff scheme is the most suitable one, 
where a three-momentum cutoff on the loop momenta is imposed on all divergent integrals. Consequently, the two- 
flavor NJL model in the chiral limit has two parameters: the coupling constant G and the cutoff A. The ultraviolet 
cutoff A determines the range of validity of the low-energy effective theory. 

We follow Hatsuda et al. [Sj , who fit the values of the four-fermi coupling G and of the noncovariant three- momentum 
cutoff A to reproduce the phenomenological parameters corresponding to the pion decay constant f.^ = 93 MeV and 
the quark condensate density per flavor (uu) = (dd) = (—250 MeV)'^ in vacuum These authors found 1231 



A = 631 MeV, G^A^ = 4.04. 



(2.3) 
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We note that at zero temperature the presence of a noncovariant three-momentum cutoff breaks Lorentz invariance 
expHcitly, hence care must be taken in implementing the cutoff in order to preserve the consequences of the chiral 
symmetry breaking in the (Lorentz invariant) full theory, viz, QCD. We will address this point in detail later in our 
discussion below. 

Following the standard Hubbard- Stratonovich procedure by introducing the composite scalar and pseudoscalar 
fields 



(2.4) 



respectively, we can rewrite £njl as 



(2.5) 



where tt = (tt^, tt^, tt"^). We note that a and tt at tree level are not dynamical fields as there are no corresponding 
kinetic terms in the Lagrangian density 1)2. 5|) . 

The focus of our study is to compare the static and dynamical critical phenomena near the critical point obtained 
from the effective field theory of scalar and pseudoscalar mesons (pions) after integrating out the quarks, and the 
0(4) linear sigma model which has been argued to describe the universality class for two- flavor QCD P, The 
Ginzburg-Landau (GL) description of the linear sigma model (LSM) near the critical point is based on the Lagrangian 
density UUE^ 



^LSM 



2 ^ 



2 



A 



B 



(2.6) 



where the coefficients A, _B > arc non-universal and vary slowly near Tc- The term emerges from the one-loop 
tadpole contribution to the self-energies which is the leading term at high temperatures and is local [23 . [25l | . 

For T < Tc we take the direction of spontaneous symmetry breaking to coincide with the scalar component a and 
split the scalar field into a mean field ctq and fluctuations S, while for simplicity we assume that there is no mean 
field along the tt direction. We write a = ao + 6. Thus the linear sigma model or Ginzburg-Landau theory in the 
approximation of mean field plus quadratic fluctuations around mean field is described by the Lagrangian density 



^LSM 

V{ao) 



V{ao) 
2 ^ 



dao 2 



cPy(ao)^2' 
dal 



f dV{ao) ^^ 
(To doQ 



(2.7) 



When the order parameter ctq is at its thermodynamic equilibrium value corresponding to dV {aa) / da^ = the tt field 
is massless and describes a triplet of Goldstone bosons. 

The goal of our study is to compare the low-energy effective field theory of meson fields obtained by integrating out 
the fermion fields in the constituent quark model, up to quadratic order in the fluctuations around the homogeneous 
mean field. This comparison will allow us to establish the validity of the universality hypothesis for both static as 
well as dynamic critical phenomena at least to lowest order in the fluctuations. 



III. STATIC GINZBURG-LANDAU EFFECTIVE THEORY 



We begin our study by obtaining the GL effective action up to Gaussian (quadratic) fluctuations around a mean field 
in the Matsubara formulation. The continuation to imaginary time that describes the theory at finite temperature 
is obtained by taking it t with < t < /3 = f/T but keeping the same Dirac gamma matrices as in Minkowski 
space-time. Furthermore, we shift the scalar condensate cr by a space-time constant mean field as 

cr(x, r) = ao + (5(x, t) , (3. f ) 

and define the effective fermion mass as 

m = Gao. (3.2) 

We have chosen the direction of symmetry breaking to coincide with the cr for simplicity, of course we could have 
chosen mean fields for the pions ttq as well. In such a case, the underlying SU(2)v isospin symmetry (for massless 



5 



quarks) entails that the mass term would be m = G\J (Tq + tTq. Choosing the a direction does not represent a loss of 
generality, by invoking the underlying isospin we can recover the general result by simply replacing CTq ^ ctq + tTq. 
The fermionic path integral can be carried out with the result for the total path integral 



Z = J VSVtv cxp I - /3V^ - J d^x dr 
-Trln(i^£; -m-G5~ iGj^r ■ tt) 



cro'5(x, r) + r) + ^^^(x, r) 



where i^E = —j'^d/dr + ^7 • V and V is the spatial volume. The fermionic path integral (partition function) 
Z[5,Tv] = cxp[Trln(i^B-m-G5-iG7V-7r)] 

DtpD^ exp J d^x j dT^{x, r) ^i^E — m ~ G5{x, r) — iG-f^T ■ 7r(x, t)] ^(x, r) 

can be written in a more familiar manner as 



(3.3) 



(3.4) 



Z[d, tt] = Zmf{ Trcxpl -G / d^x / dTtjj{S + ij^T ■ 7r)i/' 



(3.5) 

where is the r-ordering symbol and Zmf is the free-field partition function for Nc x Nf Dirac fermions of mass 
TO = Gao, namely 



Tre 



exp <^ 2pVNcN-f 



d^k 
(2^ 



-ln(l + e-'3-'=)+c^. 



(3.6) 



with ujk — Vk'^ + rn^ and (• • • )o the thermal expectation value in a free field theory of Dirac fermions with mass to, 
Nf, colors and Nf flavors. The Trln[- • •] in eq. H3.3|l can be expanded systematically in powers of G, the effective 
action in the mean field approximation plus quadratic fluctuations around the mean field is obtained by expanding 
up to 0{G^). Equivalently, the expectation value of the exponential in eq. (|3.5I) is expanded up to 0{G^) and the 
result re-exponentiated. Following these steps we are led to the finite-temperature effective action up to quadratic 
order in the fluctuations of the meson fields S and tt. Taking the spatial Fourier transform as well as the discrete 
Fourier transform in terms of Matsubara frequencies (in the interval < r < /?) 



<5(x,r) = ^^^5(k,^„)e— 



and similarly for tt*, we find 



Vn = 27rn//3, 



(3.7) 



(3.8) 



with Vn = 27rn//3 the bosonic Matsubara frequencies. 
The finite-temperature effective potential is given by 

Veff(ao) = \al - 2N,Nf 

and its derivative with respect to <tq reads 

dVoff(CTo) 



d^k [2 



(27r)3 



ln(l + e-'^'"^)+cjfc 



= (To + G{lp1p)o = (To 



1 -2NcNfG^ 



d^k 1 — 2 n{LUk) 



(27r)3 



(3.9) 



(3.10) 



where n(cjfc) = l/(e'^"''^ + 1). The equihbrium state of the system is determined by the condition dVcff{<To)/dao = 0, 
which leads to the gap equation 



(To + G{^-)p)o = (To 



1 - 2NcNfG^ 



d^k 1 - 2n{u;k) 



(27r)3 



= 0. 



(3.11) 



6 



For temperatures above the critical Tc (see below) the only available solution is uo = whereas for T < Tc there 
is a non-trivial solution which signals spontaneous chiral symmetry breaking 0| . We have chosen the direction of 
symmetry breaking along the a component of the multiplet for convenience, since in this case the tt field describes 
the multiplet of Goldstone bosons, namely pions. 
The meson self-energies are given by 

E;^(k,z.„) = -G' J d'x iyre'--^Trpl{^,T)pi{O,0)ro, (3-12) 

where the superscript c and subscript refer to the connected correlation function in the free field theory and the 
scalar and chiral densities are given by 

p(x,t) = '(/i(x,r)V'(x,T), 

P^(x,t) = ^(x,r)i7VV(x,T). (3.13) 
We discuss in detail both the effective potential and the self-energies separately below. 



A. The effective potential and its Ginzburg-Landau expansion near the critical point 



The gap equation which describes the thermodynamic equilibrium state of the system is given by 



rfVcff (gp) 

dao 



o-Q + G{'4!'ijj)o = (To 



1 ~ 2N^NfG^ 



d^k 1 - 2n{ujk) 



(27r)3 



= 0. 



(3.14) 



There is an obvious solution (Tq = corresponding to unbroken chiral symmetry, but below the critical temperature 
there is another solution with ctq 7^ which describes spontaneous chiral symmetry breaking |9j and corresponds 
to the vanishing of terms in the brackets in cq. (|3.14|) . The critical temperature corresponds to the vanishing of this 
nontrivial solution of the gap equation and is determined by the equation 



NcNfG^ 



dk k tanh - 



2r, 



27r2 



(3.15) 



where we have neglected terms that are exponentially small in the ratio A/Tc- The critical temperature is thus found 
to be given by 



V3 



NcNfG^A"^ 



A, 



(3.16) 



which for Nf — 2, Nc — 3 and the values of the parameters given by H2.3|l yields Tc ~ 150 MeV H,0. 

For T < Tf, we expand the derivative of the effective potential (|3.1U|) as a power series expansion in (Jq, this is 
achieved by the following steps: (i) In the integral in H3.10|l write the integrand as 



1 - 2n{ujk) _ 1 - 2n{k) 
Uk k 

and the integration of the first term in 13.17|l gives 



1 - 2n{uk) 1 - 2n{k) 



(3.17) 



k 1 
dk k tanh ^ - 2 



A^ 



(ii) Use the identity 



1 - 2n{uj) 



3 



{2n + l)7rT, 



(3.18) 



(3.19) 



and expand the difference in the bracket in (|3.17|) in a power series expansion in m^. The result can be integrated 
term by term. 
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Upon using the result H3.16() for the critical temperature, we find 



dm 



6 



m 



A 

2i; 



where 



^ dx 



■ tanhx. 



(3.20) 



(3.21) 



For the relevant case Nc — 3 and Nf — 2 with the values of the parameters given in (|2.3() and that for Tc given by 
eq. (j3.16(l we find F{A/2Tc) ~ 0.477. Integrating this expression we find the effective potential in the GL region (i.e., 
\T - Tc\/Tc < 1 and m/T < 1) to be given by 



Voff(TO) 



NcNf 
6 



TO' ^2 2- ™" „ / A 



2 ^ \ \2Tr 



+ 0(to^). 



(3.22) 



where the higher order terms are subdominant in the GL region near the critical temperature. In order to extract 
the momentum dependence of the fluctuation contributions to the GL effective action at finite temperature, we must 
study the self-energies. 



B. Self-energies and dispersion relations 

Anticipating our study of the real-time dynamics, in this section we summarize important aspects of the self-energies 
which will be necessary for the low-energy expansion both in the static as well as in the dynamical study. Taking the 
spatial Fourier transform, the (Matsubara) self-energies in eq. (|3.12|) are given by 

J:sik,i^n) = - dre*'^""(pk(T)p-k(0))S, (3.23) 



and similarly for E^. In what follows when no confusion occurs, we will not write explicitly the isospin indices i,j 
for to avoid cluttering of notation. The correlator in eqn. H3.23II can be written as a spectral representation 

(Pk(r)p-k(0))§ = -^^e-^^"(r.|e-^Vk(0)e--^«|TO)(m|p_k(0)|n) 

= Ve-'3^"e(^"-^'")^(n|pk(0)|TO)(m|p_k(0)|n) (3.24) 

and a similar representation for the correlator that enters in the definition of S^r h^ eq. (|3.12(l . 

The integral over the imaginary time in eq. H3.12|l can now be done straightforwardly and we find the spectral 
representation for the self-energies 

^s{k,i^n) = / dko- : , 

T,]^[k,iyn) = / dko- : , (3.25) 



where 



Ss[kM) = -^Ve-'5^"(n|pk(0)|TO)(m|p_k(0)|n)<5(fco-i?« + i?m)(e^'=«-l), 



m,n 

^2, 



S^{kM) = ■|^5]e-''^"(n|p^(0)|m)(TO|p^5_jO)|n)^(fco-i5„ + i5™)(e'"'^"-l). (3.26) 

With the purpose of studying the real-time dynamics of pion fiuctuations in linear response later and to establish 
contact between the real-time dynamics and the Matsubara formulation, let us introduce the following real-time 
correlation functions 

E>(fc,t-<') = G2(pk(i)p-k(i')>§, 

E<(fc,t-t') = G2(p_k(t')Pk(t))S, (3.27) 



8 



where pk(i) = e'^"*pk(0)e *^o*, and similarly for in terms of the chiral densities. Introducing a complete set of 
eigenstates of Hq , we find 



/oo 
duj4jk,u;)e'-^'-''\ (3.28) 
-OO 



in terms of the spectral densities 



a<(fc,w) = ^Ve-''^"(n|p_k(0)|m)(m|pk(0)|n)(5(c^ + (3.29) 

m.n 

Upon relabelling the sum indices m n in the expression for above, we find the Kubo-Martin-Schwinger (KMS) 
condition [23, Ha 

a<{k,u;) = a>{k,-u;) ^ e^^a>{k,u;). (3.30) 

Similar expressions are found for a$'{k, uj) with the scalar densities replaced by the chiral densities and the appropriate 
isospin indices 

The retarded self-energy is given by (see below) 

j:i%{k, t~t')^ ~i [E>^(fc, t - <') - S5<,(A:, t - i')] e(i - 1') = j:sAk, t - t')e(t - 1'), (3.3i) 

where the second equivalence defines Tis^-j^ik^t — t'). Using the spectral representation of the B(i — t') we find 

/oo J7 
_^g./co(*-Osret(fc^^^)^ (3.32) 
-oo 

with 

i:TUk,ko)= ' duo (3.33) 

Using the KMS condition l|3.30|l we find the imaginary part of the retarded self-energy to be given by 

Im^sAk,^) - -A^l^{k,^) - <.(fc,c^)] = ^(e''" - l)<jl^{k,u). (3.34) 
Similarly for Yis^.K{k,t — t') introduced on the right hand side in eq. (|3.31() we find 

/°° I 7 1'°° I ~ 

duje'^^'-' \<J>Ak,uj) - afAk^uj)] = - / due'^^'-' hmY.\%k,uj). (3.35) 
oo ^ J — oo 

Comparing eq. H3.34|l with a^^{k,uj) given by eq. (|3.29|) with Ss{k,ko) given by eq. 13.26|l and which enters in the 
dispersive representation of the self-energy in eq. (|3.25|l . we see that 

Ssik,ko) = linSsik,ko), 

Sl^{k,ko) = ImE;^(A:,fco), (3.36) 

Anticipating a study of relaxation in real time as an initial value problem, we introduce the Laplace transform of 
^s,TT{k,t) introduced on the right hand side in eq. H3.31|l 



^sAk.s)^ e-''^sAk.t)dt ^ — . ' . (3.37) 

Upon comparing this expression for the Laplace transform with the dispersive representation of the retarded self-energy 
H3.33(l as well as the Matsubara correlators H3.25(l . we find the following relations between the Fourier transform of 
the retarded self-energy, the Matsubara self-energy and the Laplace transform introduced above 

E^°^(fc,A:o) - Y:sAk,s^iko + e), 

T,s,Tvik,l'n) = ^S,nik,S ^ -Vn)- (3.38) 



9 



These relations will allow us to study the static aspects in the Matsubara representation and the real-time dynamics 
in an unified manner. 

The main point of this analysis and the relationship between the self-energies, is that we can obtain both the 
Matsubara as well as the real-time self-energies directly from the computation of the self-energy in the Matsubara 
representation. 

C. Calculation of the self-energies and the low-energy expansion 

Having established the relationships between the self-energies in imaginary time and in real time, the next step is 
to obtain the explicit expression of the self-energy. It is more convenient to obtain the Matsubara self-energy, once it 
is written as a dispersion relation we can use the relations given by eq. (|3.38(l and obtain the corresponding expression 
in real time. 

The most straightforward approach to obtaining the Matsubara self-energy directly in a dispersive form consists of 
writing the fermion propagator in the Matsubara representation as a dispersion relation (with zero chemical potential) 



S{p,LL!m) 



dpo- -— , w,„ = {2m + IjTrT, 



PO - lUJm 

with the spectral function for free Dirac fermions of mass m given by 



Wo 
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5{pa - Ljp) + i 



-^0 I ^-P-"^ 



5{pq +ujp) 



The one-loop self-energies are given by 

(fp 



PO - ILUm go - ii^m - iVn ' 

5 i piViPo) .5,,- p(q,9o) 



PO - lUJm 



(3.39) 



(3.40) 



(3.41) 



where q = p -I- k. The sum over the fermionic Matsubara frequencies can be done straightforwardly |24 l2fil | leading 
to the dispersive form 



(Pp 
(2^ 
n(go) - "(po) 



Y^siK'^n) = ^ y J^:^ J ^Po J d.<lQ J c?fcotr[p(p,po)p(q,go)] 

■S{ko - qo +pq) 



ko - iVn 



G 



dfco 

DO 

2 / d^P 



ImSf (fc,fco) 



(2^) 



fco — iv, 



dko 



ko - il^n 

>0 nOO poo 

3 , dpa dqo / dkotr [i^^r'^ p{p,po)ij^T^ piq, qo)] 

J —OC J —OO J ~oo 

S{ko - qo+Po) 



ImS;°t-»J(fc,fco) 
fco - iVn 



(3.42) 



where on the right hand sides of the equations above we have used the relations given by eqns. H3.37I3.38|I established 
above . 

Using the identities 



WqWp ± p • q =F = ±^[(wq ± Wp)^ - fc^ - 4m^], 

UJqUJp±p ■ q±m^ = ±-[{uJq ±UJp)'^ - k'^], 



(3.43) 
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the imaginary parts of the self-energies can be readily found to be 



lmT,f\k, ko) = -nNcNfC'^ik^ - k^ - 



-{[n{ujp) - n{ujq)][S{ka - Uq + Up) 



6{ko +Ljg- Ljp)] + [1 - n{u;p) - n{ujq)][S{kQ + ujq + ujp) - 5{ko -ujq~ t^p)]}, 

d^p 



ImE;°*'*^'(fc,fco) = ^TT5'^NcNfG^{kl-k^) 



-{[n{ujp) - n(wq)][(5(fco - t^q + iOp) 



2TT)^2LOpUJq 

-d{ko +u}q- ojp)] + [1 - ri{ujp) ~ n{ujq)][S{ka + Uq + Up) - S{ko - ujq - c^p)]}- (3.44) 

The terms proportional to [n{ujp) — n(LUq)] originate in the process of Landau damping in which a meson scatters off a 
quark in the medium, and those proportional to [1 ~n{ujp) — n{u!q)] arise from processes in which quark and antiquark 
are created via meson decay. 

It is more convenient to obtain the Laplace transform S(/c, s) for both self-energies from which the Matsubara and 
real-time ones can be obtained by simple analytic continuations. Furthermore we here focus on a long- wavelength, 
low-energy expansion in the GL region for which fc, s, to <C T ~ Tc and keep terms up to s^). After some algebra 

we obtain 



T,s{k,s) 



NcNfG^^{s^ + k^ + Afn^) 

""^ dpp^ 1 — 2n{ujp) 
— 2 1" ^ 



2 r 



dpp 
T^'^^l 



\ — 2n{ujp) \ dn(ijjp) 
Aujn 2 duj„ 



p- 



i-[l-2niu;p)] 



^'^{k,s) = S'^N.NfG'Us' + k') 



dpp^ 1 — 2 n{ujp) 



dpp^ 



dp d 
12^ dp 



1 — 2n{ujp) 1 dn{ujp) 



uj: 



Aujp 

dp d 
12^ dp 



2 duj„ 



^[l-2n(^p)] 



^ isujp is/k+p/ujp 
2kp is/k — p/tOp 



isLOp is/k + p/ujp 

1 — In — 

2kp is/k ~ p/ujp 



(3.45) 



We note that the last term in the above expressions is a surface term which arises from the introduction of the three- 
momentum cutoff A that explicit breaks Lorentz invariance even at zero temperature. If a covariant regularization 
scheme (e.g., dimensional regularization) is used, this surface term would vanish identically. The necessity of neglecting 
this surface term, a consequence of a non covariant regulator can be seen in the zero-temperature case. After the 
analytic continuation is —^ uj + ie discarding this surface term leads to the Lorentz invariant dispersion relation for 
scalar and pseudoscalar fluctuations, namely uj'^ — k'^ + Am? for the scalar {5) and cj^ = k^ for the pseudoscalar (tt). 
Hence the pions are the Goldstone bosons of the spontaneously broken chiral symmetry. Obviously the surface term 
contains an extra factor fc^ that spoils Lorentz invariance in the dispersion relation at zero temperature. Therefore 
in order to preserve the manifest Lorentz covariance, we will discard this surface term hereafter. 

The logarithmic dependence on is/k in eq. (|3.45(l is a consequence of the Landau damping contribution. Upon the 
analytic continuation is uj + ie this term leads to an imaginary part (see below) and a non-local contribution to 
the self-energy. 

D. The static effective action and Ginzburg-Landau effective theory 

The static limit corresponds to s ^ (or equivalently Vn = in the Matsubara representation). The static limit of 
the self-energies is given by 



S5(fc,0) = E5(0,0)+C(T)fc2, 
E^(fc,0) = S^(0,0) +C(T)fc2, 

where, upon neglecting the surface term in l|3.45|l . we find the coefficient C{T) to be given by 



C{T) = f ^ 



dpp^ 



2,, ,2 



1 — 2n(ujp) 1 dn[ujp) 



AlUt, 



2 dujn 



(3.46) 



(3.47) 



The long-wavelength limit of the static self-energies is given by 



^s,Tv{0,0) = hm lim 7r(fc, s). 



(3.48) 
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which are found to be given by 

£5(0,0) = -N.NfG'' 
I],(0,0) = -N.NfG'' 



^ dpp^ 1 - 2n(cjp) 



rn^ \ 1 — 2n{ujp) 2to^ dn{ujp) 

UJr, Uji dWr, 



(3.49) 



Comparing these expressions for the long-wavelcngth hmit of the static self-energies with the derivative of the 
effective potential given by eq. H3.14|l . we find the following identities 



1 + ^5(0,0) 



1 + 2,(0,0) = 



d'Veff(ao) 
dal 

1 dVoff(o-o) 



Thus the static effective action, obtained from H3.8|l by taking the - 



(3.50) 



Veff(ao) + (5(0,0 



dVcsicTo) 



+ i7r'(0,k) 



dao 

1 dVoff(cro) 



d^k 



(27r)3 



,5(0, k 



term, is given by 
d^Vofr(cro 



C(T) fc' 



5(0, -k) 



C(T) fc2 



with the effective potential in the GL region given by 



7r^(G,-k) 



4"° 



(3.51) 



(3.52) 



where the coefficients a, 6 > can be read off from eq. H3.22|l . 

Thus we see that up to a redefinition of the meson fields by renormalizing their amplitudes by the constant C{T), 
the static low-energy effective action in the GL region is exactly of the form obtained from the linear sigma model 
(|2.7|) by assuming time-independent configurations, namely setting the time derivatives in the kinetic terms in l|2.7|) 
to zero. 

The order parameter at the equilibrium minimum of the effective potential for T < Tc is given by the mean field 
solution near the critical temperature 



(r)= -(T2-T2) Cx(Te-T)^ 



(3.53) 



where we have taken the positive solution. For fluctuations around the equilibrium mean field solution in the chirally 
broken phase for T < Tc, the multiplet of pions is massless as befits Goldstonc modes, whereas the fluctuations in the 
broken symmetry direction have a mass determined by the second derivative of the effective potential at the minimum. 

This static effective action is equivalent to the expansion around a mean field up to quadratic fluctuations of the 
0(4) GL free energy 



-^gl[o-, tt] ^ d^x 



C{T) 



Vcr • Vct 



C(T) 



Vtt • Vtt 



-(T2 



TT^f +h.0.t. 



(3.54) 



where h.o.t. stands for higher order terms in <Tq,7t'^ and gradients. The equivalence between 13.54|l and H3.52|l up 
to quadratic order in the fluctuations S, tt can be immediately seen by writing a — ao + S and keeping only up to 
quadratic terms in 6,7: in the free energy H3.54II . in agreement with the discussion leading to eq. H2.7|l . Thus we see 
that long- wavelength phenomena in the static limit is in the same universality class as the three-dimensional 0(4) 
Heisenberg ferromagnet (linear sigma model) as advanced in Ref. J^, _2J . 

The local gradient expansion that allows the identification of the static effective action with the GL free energy of 
the 0(4) Heisenberg ferromagnet, is available only in the static limit. The analytic continuation of the self-energies 
to real frequencies are both non-local and complex. We now address the issue of the dynamics. 
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IV. REAL-TIME DYNAMICS 

The real-time evolution of expectation values and correlation functions is obtained from the real-time generating 
functional 

Z[J±] = TrC/(cx), -oo; J+) p U-^{oo, ~oo; J-), (4.1) 

where p is an initial density matrix and U{t, t'; J^) [U^^{t, t'; J^)] is the forward (backwards) time evolution operator 
in presence of a current J"*" (J^)- Time ordered, anti-time ordered, advanced and retarded correlation functions are 
obtained by taking suitable variational derivatives with respect to the sources J^. In the case under consideration, 
we seek to study the real-time evolution of small amplitude fluctuations around a mean field ctq, hence we take the 
initial density matrix to be thermal at temperature T < Tc in the broken symmetry phase. 

We will study the real-time evolution of small amplitude fluctuations as an initial value problem, which is defined 
by introducing external sources rjs^jr that induce an expectation value for the fluctuations around the mean field. The 
real-time generating functional l|4.1|l can be conveniently written in a path integral representation as follows 

with 

£± [ J± , ,y] = ^± (x, t)[i^-m- G(5=^ (x, t) - V • 7r± (x, i)] (x, t) + (x, t) V-^^ (x, t) 

+^±(x,<)J±(x,<) -ao^±(x,<) - i<52-±(x,t) - i7r2-±(x,i) + J±(x,i),5±(x,t) 

+J^(x,i) • 7r±(x,<) +?75(x,i)^±(x,t) + J7^(x,i) • 7r±(x,t), (4.3) 

where the superscripts ± refer to fields defined in the forward (-I-) and backward (— ) branches of the closed-time-path 
contour in the Schwinger-Keldysh formalism of nonequilibrium field theory p^ . In the Lagrangian densities above, 
the sources are auxiliary, real-time correlation functions are obtained as variational derivatives with respect to 
these and setting them to zero afterwards. In contrast, the external sources r]s,Tr above are introduced to induce an 
expectation value for the fiuctuations S, it and allow to study the real-time evolution of these expectation values as 
an initial value problem as described in detail below. 

Perturbative calculations for the meson self-energies are carried out in terms of the following real-time free fermion 
propagators: 

{rMi,''i^',t'))=z I ^Sl\t,t')e^^<--'^ , 
s;+{t,t') = S>{t,t') Q{t - t') + S<{t,t') e(t' - t), 

s^-{t,t') = sl{t,t') e{t'-t) + s<{t,t') Q{t-t') , 

S^^M^Sfit^t'), (4.4) 

where (■ • • ) denotes expectation value with respect to the initial density matrix and a,b = ±. The Wightman functions 
Sp(t, t') are expressed in terms of the spectral function as (with zero chemical potential) 

/oo 
dpop(p,Po)[l-n(po)]e-^^'''*-*'\ 
-OO 

/OO 
dpop(p,Po)n(po)e'^^'«(*-*'\ (4.5) 
-oo 

where n{po) — l/{e^P" + 1) and the spectral function p{p,po) is given by eq. H3.40|l . 

A. Linear relaxation of fluctuations 

To study the real-time dynamical evolution of the small amplitude perturbation of the a and tt fields, we consider 
preparing a system in the broken symmetry phase ctq 7^ Oj slightly perturbed away from equilibrium by adiabatically 
switching-on the external sources ry^^jr in the infinite past (see below). Once the sources are switched-off at some time. 
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say t = 0, the perturbed system relaxes towards equilibrium and the relaxation dynamics is studied in linear response. 
Since the external sources rjs^Tr induce expectation values for 6 and tt, it is convenient to decompose S and tt into an 
expectation value and a fluctuation whose expectation value vanishes (even with rfs^-^ switched-on) [2^ 

(5±(x,t) = </.(x,t)+A±(x,i), (A±(x,t))^ =0. (4.6) 

Likewise, the tt field is decomposed as 

7r''±(x,t)==xXx,t)+n*'^(x,i), (ff^±(x,t))„ =0, (4.7) 

where (x, t) is the fluctuation about the homogeneous pseudoscalar condensate which we have chosen to vanish. 
The expectation values (n^(x, i))^ and (A^(x, i))^ are obtained consistently in a loop expansion with the real-time 
generating functional in presence of the external sources rj. 

The equations of motion for the expectation values in linear response are obtained implementing the tadpole 
conditions (11='= (x, t))^ = and (A='=(x, t))^ = 0. These expectation values are obtained in a systematic loop expansion 
in terms of quark loops but only to linear order in (/)(x, t); x*(x, t) consistent with linear response. This method has 
been implemented in many different problems, and the reader is referred to reference|28i| for more details. 

The linearized equations of motion for the spatial Fourier transform of the small amplitude condensate fluctuations 
(/)(x, t) and x(Xi^) ^-re the following 

</.(k,i)+ /* dt'j:sik,t-t')^ik,t') + ^^^«(^"^ j(3)(k) ^ 7y,(k,t), 

J-oo dao 

x\Kt)+ f dt'j:^{k,t-t')x\Kt') = vUKt), (4.8) 

<y —oo 

respectively. Since the quark fields are described by a thermal initial density matrix, the real-time self-energies 
S5,7r(fc,i — t') are related to the respective retarded self-energies of the a and tt fields given by eqs. (|3.31l) - (|3.33l) and 
(|3.44l) . In this section we consider fluctuations around the equilibrium state, namely (iVcff((7o)/(icro = 0, postponing 
the discussion of fluctuations off the equilibrium minimum to a later section (see sec. IIV C|) . 

The adiabatic preparation of a system slightly perturbed away from equilibrium can be realized within the real-time 
formulation by taking the spatial Fourier transform of the external sources to be of the form p9| 

mA'^,t)^mA'^)e''Qi-t), e^0+. (4.9) 

The e-term serves to switch on the source adiabatically from t = —oo so as not to disturb the system too far from 
equilibrium in the process. If the system was in an equilibrium state at t = — oo, then the equilibrium condition 
dVeff{(To)/dao = ensures that for t < there is a solution of the equations of motion l|4.8|) of the form 



0(k,t) =0(k,O) e^* 



,t for t < 0, (4.10) 



where 0(k, 0) [x*(k, 0)] is related to ?7i(k) [77^ (k)] through the equations of motion for t < 0. That this is a consistent 
solution for i < can be seen easily from the representation H3.35|l . 

The integrals over t' in the non-local terms in the equations of motion 14.811 can be separated into t' < Q and t' > 0, 
in the interval — 00 < i' < the solution is given by eq. H4.10|l . There only remains a convolution for t' > which 
is studied by taking the Laplace transform of the equations of motion for t > 0. The equations of motion for the 
fluctuations, in terms of the Laplace variable s, are given by 

DAik,s)^ik,s) = i [5,-1(^,5) -5,-1(^,0)] 0(k,O), 

D-\k,s)x'{k,s) = -^[D-\k,s)~D-\k,0)]x\KO), (4.11) 

where 0(k, s) [x*(k, s)] is the Laplace transform of (f>{'k,t) [x(k, t)]. In the above equations, D'^^(k,s) is the inverse 
retarded propagator of the corresponding field, namely 

5,-i(fc,5) = l + E5,,(fc,s), (4.12) 

where Yis,-n{kTs) is the Laplace transform of the self-energy given by equation (|3.37|l . In the long- wavelength, low- 
energy limit, these Laplace transforms are given by equation (|3.45|) . The self-energy I]5.Tr(fc, lo) and inverse propagator 
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Dg ^{k, uj) in term of the real frequency w are defined as the boundary values of the analytic functions E5^^(fc, s) and 

Dj^{k, s), respectively, through the analytic continuation s —>■ iuj + e. 

To simplify the notation, we will henceforth drop the isospin indices of the pion field by invoking the remaining 
SU(2)v isospin symmetry. The real-time evolution of the small amplitude fluctuations ^(k, t) and x(k, t) for t > is 
now in the form of an initial value problem with the respective initial value specified at <: = given by 0(k, 0) and 
x(k, 0). The solution of this initial value problem can be obtained from the inverse Laplace transform, e.g., [similar 
expression for (/)(k, t)] 



ds 



(4.13) 



where x(k, s) is the solution of (|4.11|l and the Bromwich contour B runs parallel to the imaginary axis in the complex 
s-plane to the right of all the singularities (isolated poles and branch cuts) of x(k, s) . We note that there is no 
isolated pole in (/)(k, s) and x(k, s) at s = since the corresponding residue vanishes. 

We now focus our attention to the real-time dynamics of small amplitude fluctuations near the equilibrium minimum 
determined by dVcs{'^o)/d(Jo = 0. It proves convenient to introduce the following dimensionless variables 



k = 



s 



the inverse propagators of the a and tt fields now read 



T' 



A 



(4.14) 



D-'(k,s) 
D-\k,s) 

where we have used eq. (|3.14ll and introduced 



g{s^ ^k^ + Ax^)H{iS_,x), 



(4.15) 



''{z,x) 



tanh 



1 



e{z, x) 



1 - 



i^e{z,x) ^ i£^ + z/e{z,x) 



2z 



i^ — z/e{z, x) 



e{z, x) = + ; 



2cosh2£(i£) 

(4.16) 

The (quasi)particle excitations of the scalar meson are massive and as discussed in Ref. 8] they correspond to 
a broad resonance arising from the decay of this meson into quark-antiquark pairs even at lowest (one-loop) order. 
Physically the scalar excitations are not stable. In what follows we will restrict our attention to the dynamics of pions 
which are the true low-energy degrees of freedom in QCD. 

While the one loop self-energies for the scalar and pseudoscalars are fairly well known "^j*^ and to this one loop order, 
or alternatively leading order in the large N^Nf corresponds to the random phase approximation (summation of ring 
diagrams for the propagators), unlike previous studies, we focus on the real time dynamics for which the formulation 
described above is necessary. The analysis of the static case, while in agreement with previous studiesj^ allows 
to study on the same footing and under the same approximations both the static as well as dynamical aspects thus 
allowing an unambiguous comparison. In particular the explicit form of the static effective action obtained above 
will allow us to compare the predictions for the growth rate of spinodally unstable long-wavelength fluctuations as 
obtained from the static effective action and those obtained from a detailed analysis of the dynamic response. 



B. Pion dynamics 

Fluctuations around equilibrium 



We begin our study by focusing on the real-time (quadratic) fluctuations around the equilibrium broken symmetry 
state for T < Tc, which is determined by the condition dVcf{{(To)/daQ = with ctq 7^ 0. In this case, the real-time 
dynamics of pions is determined by the Laplace transform 



X(k, s) 



H{0) 



X(k,0). 



(4.17) 



After the analytic continuation s 



VjJ ■ 



ii+e)H{itx)_ 

e and introducing the dimensionless variable a = iv/k, we find the real 
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FIG. 1: The function H-Ei{a,x) vs a for different values of x. 



and imaginary parts of the analytically continued function H{a + ie) given by 



Hnia, x) 



47r2 



dzz 



2 r 



tanh 



e(z,x) 



e(z, x) 



2cosh2l(^ 



Hi{a,x) = 



2 

max 



dzz 



respectively, 



Q:£{z, x) 
2i 



In 



a + z/e{z, x) 



a — z/e{z, x) 



,i„ e{z,x) cosh^ £(^£) 



(4.18) 



"max = A/ VA2+^ ; z^in = a;a/ Vl - q;2 (4.19) 
The real part can be written in the following manner which separates the terms that are most sensitive to the cutoff 



B{a, x) 



1 



1 



47r2 



dz 



tanh + z^) 



1 



dzz^ 



e2 2cosh2^ 



,3/2 



ae(z,x] 



-e(z-i) 



i3(a, a;) 



2z 



In 



a + z/ e{z, x) 



a ~ zj e(z^ x) 



(4.20) 
(4.21) 



The second term in eqn. (|4.2U|) is rather insensitive to A and in this expression A can be taken to infinity with an 
exponentially small error. The form of i?(a, a;) above displays explicitly the dependence on A. 

The real and imaginary parts given in (|4.18l) are plotted in Figs. ^ and |21 respectively, for Nc — 3, TV/ — 2, the 
values of G, A given by eqn. (|2.3|l and different values of x in the Landau-Ginzburg region x 1. 

The real part i?R(ck,x) increases monotonically for < a < amax then decreases monotonically to its asymptotic 
value for a 3> 1. Furthermore, HB,{a,x) features a sharp cusp near a < amax as x — ^ 0. This cusp implies infrared 
divergence in the self-energies at the critical point a; = and, as will be seen below, has a importance consequence 
on the relaxation of critical fluctuations. The imaginary part Hi{a,x), which is completely determined by Landau 
damping only has support in the interval —amax < a < amax- 

The analytic properties of the quark loop diagram at finite temperature had also been studied in ref . |3(| . 

The analytic continuation of the propagator to real frequency features the following singularities: 

• Isolated poles. Isolated real poles describe stable quasiparticle excitations, whereas complex poles (on the 
second sheet) correspond to resonances with finite width. From H4.11|l and H4.15(l . isolated poles of the Laplace 
transform for the pion field are determined by the equations: 



(1 - a^)HR{a, x) = 0, (1 - a^)Hi{a, x) = 0. 



(4.22) 
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FIG. 2: The function Hi{a,x) vs a for different values of x. 



Since Hi{a,x) vanishes at a = ±1, there are isolated real poles at a = ±1 provided that HB,{a,x) remains 
finite there. This indeed is the case below the critical point when x > 0. Specifically, 7?p{(l) is a monotonically 
decreasing function for a; > and its asymptotic behavior for x ^ 1 is found to be given by 



In 



(4.23) 



Thus for X > pions are sharp quasiparticles in this approximation that propagate with a dispersion relation 
^Trik) = ±fc, namely with unit group velocity. 

• Branch cuts. Branch cuts represent multiparticle states. In the long- wavelength, low-frequency limit there 
is only a branch cut on the real frequency axis in the region — amax^ < uj < amax^, corresponding to Landau 
damping processes in which the pion scatters off a quark in the medium. 

Having identified the singularities of the Laplace transform x(k, s), we can perform the inverse Laplace transform 
(I4.13|) and obtain the real-time evolution 



X(k,t) 
X(k,0) 



da 



-p{a) cos ar, 



(4.24) 



where r = kt, p{a) is the spectral function for the pion. In what follows we set x(k, 0) = 1 for notational simplicity. 
The spectral function p{a) receives contributions from the quasiparticle poles and the Landau damping cut: 

p{a) = Zsgn(Q;) S{1 - a^) + Pcut(a), (4.25) 

where sgn(Q;) is the sign function, Z is the residue of the pion poles at a = ±1, and Pcut(a) denotes the cut contribution 



Hn{0,x) 
^ = ^^^T r, Pcut(a) = 



HniO,x)Hi{a,x) 



Tr{l - a-^)[Hl{a,x) + Hf{a,x)y 



Hk{1,x)' 

Therefore the real-time solution is given by 

X(k, t) = Xpoie(k, t) + Xcut(k, t), 

where for x(k, 0) = 1, we find 

Xpoio(k, i) = Zcosr, 

Xcut(k, i) = 2 /" — pcut (a) cos ar, r = kt. 
Jo Oi 



(4.26) 



(4.27) 



(4.28) 
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FIG. 3: The residue of the pion pole Z vs x. 



Setting < = in H4.28|l [with x(k, 0) = 1], we obtain the following sum rule for the spectral function 



^0 



Z + 2 



— Pcut(a) = 



(4.29) 



a 



which we have verified numerically for a wide range of x. 

The quasiparticle poles give rise to an undamped oscillation with an amplitude determined by the residue Z . 
Figure 13 displays the residue of the quasiparticle pole Z as a function of x. It shows clearly that away from the 
GL region in the low-temperature limit (a; ^ 1) the spectral function is dominated by the isolated pion pole since 
Z K. \. Thus Landau damping is negligible and pions are stable excitations. This is because Landau damping is 
solely a medium effect arising from scatterings of pions with quarks in the heat bath, hence its contribution at lower 
temperature is less important. In the GL region (a; <C 1) the continuum contribution from the Landau damping 
cut dominates the spectral function, since in this region Z <^ 1. Figure |21 also shows a crossover at a; ^ 1 at which 
Z ~ 1/2, hence the isolated pole and the continuum contributions are of equal importance. The Landau damping 
contribution to the spectral function jOcut(<^) is displayed in Fig. 01 For a; ^ 1, it reveals clearly a sharp peak near 
the threshold, indicating the emergence of an infrared divergence at the critical point. 

Equation H4.23|l and the fact that H^ {0) is finite and non- vanishing at a; = entail that near the critical temperature, 
namely x <ti 1, the residue of the pion pole behaves as 



where we have used that x = m{T)/T and m{T) = Gaeq{T) cx (T^ - T)^/^ by eqn. (|X551) . 

Therefore the residue at the isolated pion pole vanishes as T — > from below. As the critical point is approached, 
the weight of the isolated pole to the spectral function becomes negligible and exactly at the critical point the residue 
vanishes and the isolated pole merges with the continuum contribution. Thus it is clear that as the critical region is 
approached, the edge of the Landau damping cut begins to dominate the spectral function. 

The Landau damping contribution to the real-time evolution of a long- wavelength pion fiuctuation Xcut (k, r) is 
plotted in Fig. El The value at r = equals 1 — Z by the sum rule, thus it is clear from this figure that for small x, 
namely T < Tc, the continuum contribution dominates over that of the isolated pole. This figure also reveals several 
time scales: a rapid damping as a function of the scaled variable t = kt towards an oscillating amplitude, which is 
of the same order as the isolated pole contribution Z <C 1 and oscillates with a frequency a ~ 1 and slowly damped 
out on longer time scale . The frequency of these oscillations is clearly determined by the sharp rise of the Landau 
damping contribution near a = 1. 

Although for T < Tc the pion is a stable quasiparticle, its oscillating amplitude Z vanishes as the critical temperature 
is approached thus near the critical point the initial amplitude is completely damped out by Landau damping. As 
T — > Tc from below the amplitude of long-wavelength pion fluctuations is damped by Landau damping on a scale r, 
which is numerically insensitive to x for a; < 0.2 and is approximately ~ 5 for the values of coupling and cutoff 
given by eq. (|2.3|) . In terms of the physical time, we see that the amplitude of long- wavelength pion fluctuations 
relaxes on a time scale trci(fc) = T^/k^ where r* ~ 5 is numerically found to be rather insensitive to x. Thus relaxation 
of long-wavelength pion fluctuations near the critical point are critically slow. 
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FIG. 5: Xcut(k, r) vs r in the GL region. 



2. Dynamics at the critical point 

As remarked above, the self-energies feature a logarithmic threshold divergence at the critical point. Indeed, for 
X — > the real and imaginary parts of H{a + ie) can be approximated, respectively, as 



dz 







z 




a4 
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\ u ^1 




(l-f.n 




tanh — 


2 ^ 
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/ 2_ 



^6(1 



2N , A 
Y tanh—. 



(4.31) 



Whereas the analytically continued inverse pion propagator vanishes as — > 1 (with a logarithmic divergent slope 
for the real part), the values a = ±1 (w = ±A:) are longer identified with quasiparticle poles but the endpoints of 
the logarithmic branch cut. This threshold singularity is reminiscent to those in QED and in equilibrium critical 
phenomena ,21.]. In the former case the self-energy of a charged particle has a logarithmic infrared divergence near the 
threshold arising from the emission of soft gauge bosons. In the latter case logarithmic infrared divergences appear 
in correlation functions near the critical point at the upper critical dimension. The emergence of infrared divergences 
in both cases indicates the breakdown of perturbation theory and calls for a nonperturbative resummation of the 
infrared divergences [22. 

To display the logarithmic threshold divergence in an explicit manner, let us focus on the Laplace transform of the 
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inverse pion propagator [see (|4.15() ]. which near the threshold at the critical point can be approximated as 

D-\k,s) = gCP{l+e) 

gCP{l+e){l-^H-{^ + e)] + ---}, (4.32) 

where 

C = Fr(0), a = ^ tanh | « 0.683, (4.33) 

and the dots denotes terms that are formally of order A but remain finite in the limit — )■ —1. Such logarithmic 
threshold divergence at the critical point is akin to that found in Ref. [23|. In this reference a finite-temperature 
version of the renormalization group was introduced to provide a resummation of these threshold divergences which 
leads to an exponentiation and anomalous dimensions near threshold. 

The logarithmic threshold divergence featured by the pion propagator H4.32|l suggests the buildup of an anomalous 
dimension. The results of Ref. 0| would suggest that the leading logarithmic enhancement of the pion self-energy at 
threshold eq. H4.32II could exponentiate leading to a resummed inverse pion propagator near the pion mass shell 

D-\k,s)c,-gCP[-{l+e)]''^ for ^ -1, (4.34) 
which, after analytic continuation iS^ a + ie, gives 

D-i(fc,cj) ~ -(/CP |l-a2|i-^{cos[A7re(l-a2)] 

+isgn(a) sinA7re(l - a^)} for 1. (4.35) 

Such an exponentiation if correct would reveal a new dynamical exponent (anomalous dimension) A, which however 
is numerically large to be justified by a perturbative approach. 

Obviously at this stage we can only speculate that such an exponentiation actually emerges based on the study 
of Ref. [2l[. A full analysis of whether a finite-temperature renormalization group approach can be systematically 
implemented in this theory, remains at this stage an open question that deserves further study and which is clearly 
beyond the scope of this article. 
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3. Non-locality of the real-time effective action 

An important aspect of the pion self-energy is that it is non-analytic for fixed /c as — > 0. In particular this 
non-analyticity of the function H{a) entails that the limits fc — > and — )■ rfo not commute. Furthermore, this 
branch cut prevents a well-defined expansion in powers of A: or cj, which in turn translates in that the effective action 
is non-local in space and time and cannot be expressed as a gradient expansion in terms of space and time derivatives 
of the fields. This non-locality is similar to that in the hard thermal loop effective action in gauge theories |2^ |^ 
which also originates in Landau damping. 

To highlight this point more clearly, we focus on the limits A;, a; ^ of the retarded inverse pion propagator [see 

eqs. (gunn 

D-^\{k,uj) = 5^\fc,s = iuj + e)^ ~{uj^ - P) [H^(a,x) + iHi{a,x)] . (4.36) 
We find the following limits 

D;^,{k,Q) = eH^{Q), D-J(0,c^) = -c^2^R,(oo), (4.37) 

with H^ifQ) and H^{oo) slowly varying functions of the temperature near Tc- Their values at T — Tc are given by 

i/R(0,r = r,) = 0.015, i7R(oo,r-Te) = 0.04. (4.38) 

The limits (|4.37|l would seem to suggest that there are stable long-wavelength pion fluctuations that propagate with 
a dispersion relation 



u;^^v'k', „2 ^ ^r(0) ^ 375_ (4 39) 

Hk[oo) 

However, such conclusion is obviously incorrect, for T < Tc there are stable, long-wavelength pion fluctuations that 
propagate with unit group velocity since the isolated pion poles correspond to = 1 [see eq. (14.22(1 ]. Landau damping 
thus prevents a local description of the low-energy effective action as an expansion in spatiotemporal gradients. This 
point will be important below when we discuss our results in light of previous studies. 
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C. Away from equilibrium: spinodal instabilities 



During a rapid phase transition, early stages of the nonequilibrium dynamics is dominated by spinodal instabilities, 
i.e., long- wavelength fluctuations are unstable and their amplitude grows exponentially during the early stages [25| . 
These instabilities result in the emergence of correlated domains, which in the case under consideration are coherent 
domains of pions which have been identified as disoriented chiral condensates . Spinodal instabilities emerge when 
the order parameter is not at the minimum of the free energy for T < T^, but is still "rolling" towards it. 

In the case under consideration the equilibrium minimum in the GL region is determined by eq. H3.53|l . Hence 
below the critical temperature and in the GL region, we have 

^dV^ = b [4 - . »A^(r), (4.40) 

(To CtCTo 

with A2(T) > in the spinodal region |cro| < crcq(T). The inverse pion propagator (|4.15|) in the GL region away from 
equilibrium can now be written as 

D-\k, s) = g{s^ + ~k^)H{ii, x) - A2(T), (4.41) 

where we have restored the term in (|4.40() which is absent in equilibrium. It is clear from the inverse Laplace 
transform expression that yields the real-time evolution H4.13(l that real poles in the Laplace variable s correspond to 
exponentially growing (or decaying) solutions. If the real part of the pole is positive, the real-time evolution describes 
an instability with a consequent exponential growth of small amplitude fluctuations. 

For spatially homogeneous pion fluctuations (fc = 0) in the spinodal region of the phase diagram with |cro| < (Teq(T) 
for which /\^(T) > 0, the inverse pion propagator becomes 



D~\0,s)^gs^I{T)~A^{T), (4.42) 



where 



1 dzz^ e , , 1 /2 A e'^ , , , 

^ 4^ / ^ 2 ' ^(^^^ ^ 4^ ( ^ ■ (^-^^^ 



with 7 the Euler-Mascheroni constant. The pion propagator has real poles at 



which describes one growing and one decaying solution. The growing solution has growth rate given by (for T < Tc) 

r,(fc = 0,r) = (X [al^{T) ^ alY . (4.45) 

The classical spinodal corresponds to the region of order parameter in the phase diagram for which the system 
is thermodynamically unstable to long- wavelength perturbations 25] . This corresponds to the values of the order 
parameter for which the zero momentum growth rate r,(fc — Q,T) ^ 0. Therefore we see that the classical spinodal 
line corresponds to |cro| < crcq(T'). This is the same as the classical spinodal in the mean field approximation for the 
GL effective theory or the linear sigma model. This can be seen from the Lagrangian density (|2.7() . Since along the 
classical spinodal (I/ctq) dV {<tq) / da^ < 0, hence the mass squared term for the pions is negative and small amplitude 
pion fluctuations are unstable and grow exponentially 25]. 

While obtaining the growth rate for arbitrary k is an involved numerical task which is not very illuminating, we 
can learn much by obtaining T^,{k,T) in the long-wavelength limit k'^/A'^{T) <C 1. This is achieved by expanding 
the function H{is/k) for large values of the ratio s/k and solving for the position of the (real) pole in the variable 
s. This expansion is valid since the function H is analytic everywhere away from the Landau damping cut along the 
imaginary axis in the complex s-plane. After straightforward algebra we find 



r*(/c,T) 

where 



/\^(T^ k'^ / 

""^'^ [l^J[T)]e+gI[T)J{T)-!^^+o( 



gI{T) ' ' " ' a V / V ^^2(y) VA4(r) 



(4.46) 



247r2yn cosh^ ^ ' Utt^ 



JiT)^^J -^^:TT ; ^(T^c)-T^. (4.47) 



2 



21 



While the expression H4.46|l can be expanded in a power series in fc^ in the long- wavelength limit, we have purposely 
kept the powers of k inside the square root to compare with the growth rate predicted by the linear sigma model l|2.7|) 



(4.48) 

The corrections to the coefficient of fc^ as well as the higher powers of k in the growth rate ()4.4tif) all originate in the 
non-local Landau damping contribution. 

V. CONCLUSIONS AND SUMMARY OF MAIN RESULTS 

In this article we have studied the low-energy effective action for meson fields in the two-flavor NJL constituent 
quark model. The main goal and focus are to study static and dynamical phenomena near the critical point for the 
chiral phase transition and explore whether the low-energy effective action describes the same universality class as 
the 0(4) linear sigma model Q or the nonlinear sigma model [2^ . 

We have obtained the low-energy effective action by integrating out the quark fields to lowest order in the loop 
expansion up to quadratic order in the fluctuations around a mean field configuration. In this constituent quark 
model, the meson degrees of freedom emerge as scalar and pseudoscalar bound states of quark and antiquarks and 
their space-time dynamics is completely determined by the underlying quark model. Both the static as well as the 
dynamical aspects where studied in the low-energy, long-wavelength limit and in the GL region near the critical point 
in the broken symmetry phase. The evaluation of the effective potential (free energy for a mean field configuration) as 
well as the self-energies for the scalar and pseudoscalar mesons allow us to obtain the effective action from which we can 
extract the static limit. The scalar and pseudoscalar self-energies are obtained to lowest order in the loop expansion 
upon integrating out the quark degrees of freedom. The imaginary parts of the self-energies in the low-energy, long- 
wavelength limits are completely determined by Landau damping resulting from quark (antiquark) scattering in the 
medium. The integration of the quark degrees of freedom in the real-time effective action allows us to obtain the 
equations of motion for scalar and pseudoscalar (pion) meson fluctuations around the mean field and to study the 
real-time evolution of small amplitude fluctuations as an initial value problem. The study of the real-time dynamics 
near the critical point allows us to explore dynamical critical phenomena. 

Our conclusions and main results are the following: 

(i) Static critical phenomena. We have established that the low-energy effective action in the static limit up to 
quadratic order in the fluctuations around the mean field, which is obtained by integrating out the quark fields 
up to one loop, coincides with the GL free energy including the spatial gradients in the long-wavelength limit. 
Thus the static limit of the NJL constituent quark model with two flavors is in the same universality class as 
the 0(4) linear sigma model (Heisenberg ferromagnet) in agreement with the conjecture in Ref. 0- 

(ii) Dynamical critical phenomena. Dynamical phenomena of pions near the critical point is studied by focusing 
on the equations of motion for small amplitude pseudoscalar fluctuations around a mean field configuration. 
The equations of motion are solved by Laplace transform as an initial value problem. For fluctuations around 
the equilibrium minimum of the effective potential below the critical temperature, we find that there are stable 
long-wavelength pion excitations corresponding to poles away from the continuum in the pion propagators. The 
dispersion relation of the stable pion excitations is uJT,{k) = k and the residue at the pion pole vanishes as 
Z (X 1/1 ln(l — T /Tc)\ as T ^ Tc from below. The amplitude of the pion fluctuations is Landau damped on 
a time scale trci(fc) — T^^jk with r, a slowly varying function of temperature near T^. Thus near the critical 
temperature relaxation of long-wavelength pion fluctuations are critically slow. At the critical temperature, the 
isolated pion pole merges with the Landau damping continuum resulting in a logarithmic divergence on the 
pion mass shell similar to the logarithmic enhancement found in Ref. . Based on the renormalization group 
arguments of Ref. we conjecture that such a logarithmic divergence on the pion mass shell at the critical 
temperature is the harbinger of a dynamical anomalous dimension A, which however turns out to be fairly large 
in this model, A ~ 0.683. 

We have studied the dynamics of fluctuations away from the equilibrium minima of the effective potential to 
assess the growth rate of spinodal fluctuations. We find that the classical spinodal line coincides with that of 
the mean field description of the GL linear sigma model. However, the growth rate for long-wavelength spinodal 
fluctuations is different from those obtained from the GL effective theory. The Landau damping contribution 
to the self-energy leads to a richer wavelength dependence of the growth rate of small amplitude spinodal 
fluctuations. 
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Thus we conclude that whereas static critical phenomena is described by the universality class of the 0(4) Heisenberg 
ferromagnet (linear sigma model) near the critical temperature, dynamical critical phenomena is not in the same 
universality class. We also conclude that below the critical temperature there are stable pion excitations with a 
dispersion relation w^Cfc) = k, namely with unit group velocity. This result is in disagreement with those found in 
Ref. ||22|. We argue that the source of disagreement is likely to be found in the derivative expansion of the effective 
action as proposed in Ref. 22]. We argue that such derivative expansion leading to a local effective action is not 
reliable because the Landau damping contribution to the pion self-energy introduces a cut in the complex frequency 
plane along the real axis, which at the critical point is in the region —k < to < k. This cut results in that the 
real-time effective action is non-local and cannot be expanded in spatio-temporal derivatives. The results of our one 
loop calculation must be compared with the results of the one-loop calculation in the linear sigma model in the first 
article of ref.|22 where the limits fc ^ 0, w = and the limit a; ^ 0, fc = are studied. The Landau damping cut for 
—k<uj<k where the pion propagator is non-analytic makes these two limits not commutative. We suspect that this 
non-analyticity is at the heart of the discrepancy between our results and those of ref.|22, at least at the one-loop 
level which is common to both studies. 

Such a non-analytic structure below the light cone is akin to that found in the hard thermal loop limit of gauge 
theories which also arises precisely from Landau damping |i26i. i3Lj . At the critical point we find a logarithmic divergence 
on the pion mass shell , ag ain originating in the Landau damping processes. These logarithmic divergences are akin 
to those found in Ref. 21] and suggest the buildup of (in this case a large) dynamical anomalous dimension. 

Caveats: beyond lovifest order. Our results are based on integrating out the quark degrees of freedom at lowest 
order in the loop expansion. Given that the effective coupling is not small the validity of such an expansion can 
be called into question. Such expansion becomes exact in the \a,rge-NfNc limit, however, taking NfNc to be large 
for fixed coupling one can see from eq. (I3.16f) that Tc ~ A and the regime of validity of the constituent quark 
model as a description of the critical properties must be re-assessed. In the linear (or nonlinear) sigma model effective 
field theory, the self-energy corrections arise from interaction terms between scalar and pseudoscalar fluctuations. 
For example some of the results of Ref. [22] are based on a one-loop contribution to the pion self-energy from the 
S-TT-TT vertex below Tc- One could obtain the effective J-tt-tt vertex from the constituent quark model also in one-loop 
approximation by going beyond quadratic order in the fluctuations. Such vertex necessarily must be proportional to m 
by chiral symmetry. A local limit of such vertex can be obtained by taking first the static and then the long-wavelength 
limit. However, using such local vertex to calculate real-time dynamics would not be consistent and the non-local form 
factor of the vertex must be taken into account to calculate the corrections to the pion self-energies. Such calculation 
is clearly much more complicated than that in the local effective field theory because of the momentum and frequency 
dependence of the vertex are likely to change the results of Refs. I2lll2^ . 

At the order considered, below the critical temperature the pion pole is isolated, this is a consequence that to this 
leading order there are no scattering contributions and therefore the pion does not acquire a collisional width. Of 
course this situation will change beyond this order, in particular corrections to next to leading order in the large 
NcNf limit will lead to scattering corrections and the pion will acquire a collisional width, whose behavior near the 
critical point must be studied. The collisional width of the pion is an important ingredient in the manifestation of the 
axial anomaly at finite temperature js^] whose understanding thus requires to go beyond one-loop order, which is the 
leading order in the large NcNf limit. Furthermore, as we mentioned in the introduction, a well known shortcoming 
of the NJL model is its lack of confinement [^Q, which is manifest in that free quarks are present in the thermal bath 
below the critical temperature for chiral symmetry. Clearly there are no models, short of full QCD that can overcome 
this problem. 

The logarithmic threshold enhancement near the critical point suggests that a resummation scheme such as a 
renormalization group program as proposed in ref. j2ll] must be implemented to approach the critical region in a 
reliable manner. 

The study of the long-wavelength small frequency limit requires a resummation of self-energy and vertex diagrams, 
just like in the case of transport phenomena 33, 34] as a consequence of the emergence of "pinch singularities" js^ 
which are manifest as secular terms in the perturbative expansion in real time|34|. In ref. |34j it was shown that such 
resummation can be implemented via a dynamical renormalization group. It is clearly interesting and important to 
study the long-wavelength and small frequency limit near the critical point implementing these resummation schemes, 
perhaps supplemented by a renormalization of the scattering amplitudes near the critical point as suggested in ref. [2l] . 
Clearly such program is well beyond the one loop approximation which despite the fact of being the leading order in 
the large NcNr limit, does not include self-energy and vertex corrections which are necessary for a full resummation 



The implementation of such program is beyond the scope of this study. 

Thus our study suggests that the static critical phenomena is well described by a GL low-energy effective theory, but 
that the dynamical aspects are richer, with novel dynamical critical phenomena whose universality must be studied 
in detail. Our study at the one-loop level should be taken as indicative of novel dynamical phenomena emerging near 
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the critical point and calls for resummation program perhaps by implementing a renormalization group approach akin 
to that of Ref. (21| but directly in the constituent quark model for a deeper understanding of dynamical scaling and 
anomalous dimensions. 
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